A unified approach to split absorbing boundary conditions 
for nonlinear Schrodinger equations 

Jiwei ZhangJ*| Zhenli XuJ3 and Xiaonan Wil] 
Abstract 

An efficient method is proposed for numerical solutions of nonlinear Schrodinger equations in 
an unbounded domain. Through approximating the kinetic energy term by a one-way equation 
and uniting it with the potential energy equation, absorbing boundary conditions are designed to 
truncate the unbounded domain, which are in nonlinear form and can perfectly absorb the waves 
outgoing from the truncated domain. We examine the stability of the induced initial boundary 
value problems defined on the computational domain with the boundary conditions by a normal 
mode analysis. Numerical examples are given to illustrate the stable and tractable advantages of 
the method. 
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I. INTRODUCTION 



In this paper, we consider the construction of absorbing boundary conditions (ABCs) for 
the nonlinear Schrodinger equation (NLS), 



ih 



dip(x, t) 



h 2 d 2 



v(x) + /(H 2 )U, 



(i) 



2m dx 2 

where m is the atomic mass, h represents the Planck constant, and V(x) is the potential 
function. This equation appears in many different applications Q, sucn as ^ ne gravity 

waves on deep water in fluid dynamics, the pulse propagations in optics fibers, and the Bose- 



Einstein condensation. For the cubic nonlinear term f(\ip\ 2 )ijj 
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ifr, Eq. (pQ) reduced 



12J and references 



to Gross-Pitaevskii equation, which be numerically studied by 
therein. For a quintic nonlinearity, one can see [l7]. 

A current challenge for the numerical solutions of this kind of problems is due to the 
unboundedness of the physical domain. The common practice to overcome this difficulty is 
to limit the computational domain and solve a reduced problem. To make the truncated 
problems complete, boundary conditions should be added to the reduced problem. ABCs 
applied at prescribed artificial boundaries have been widely studied in recent decades, such 



as 



23] for hyperbolic wave equations, |2l|, [381] for el 



iptic equation s. |16l . 1201 ] for parabolic 



equations; see also the reviews by Givoli[9(, Tsynkov[3j|, Hagstrom|l5| and Han|l9j. The 
purpose of designing ABCs is to annihilate all the incident waves so that minor reflected 
waves may propagate into the computational domain. ABCs can be distinguished into 
two categories: the global ABCs and local ABCs. The global ABCs usually lead to the 
well-approximated and well-posed truncated problems, but the implementation cost for the 
global ABCs is expensive. Fast evaluation methods 24| should be applied to treat the global 
ABCs. On the other hand, local ABCs are computationally efficient, but the accuracy and 
stability are the main concerns. There is another way to design ABCs based on the media 
of the material, called the perfectly matched layer (PML) methods {3], which have been 
applied to many complicated wave propagation problems. 

For the nonlinear equations, however, it is difficult to find a suitable boundary condition 
on the artificial boundary in general, except for some special cases such as the problem can 



be linearized [18j, or a Burgers type questions [22 



34 



37| which can be corresponded with 



the linear parabolic equations by the Cole-Hopf transformation. For the NLS equation on 



the unbounded domain, Zheng 



321 ] obtained the transparent boundary condition using the 



inverse scattering transform approach for the cubic nonlinear Schrodinger equation in one 
dimension. Antonie et al. [2| also studied the one-dimensional cubic nonlinear Schrodinger 
equation and constructed several nonlinear integro-differential artificial boundary conditions. 



In 



29 



331 ]. Szeftel designed absorbing boundary conditions for one-dimensional nonlinear 
wave equation by the potential and the paralinear strategies. Soffer and Stucchio]28] pre- 
sented a phase space filter method to obtain absorbing boundary conditions. The PML 



39 



was also applied to handling the nonlinear Schrodinger equations. The recent pa- 



pers [35|, |36|, split local absorbing boundary (SLAB) conditions have been developed through 
a time-splitting procedure for one-dimensional nonlinear Schrodinger equations. The local 
absorbing boundary conditions were imposed on the split linear subproblem and yielded 
a full scheme by coupling the discretizations for the interior equation and boundary sub- 
problems. In addition, an adaptive method was developed to determine the parameter ko 
involved in the local boundary conditions. 

In this paper, we present an efficient implementation of the nonlinear absorbing boundary 



conditions for the NLS equation by making use of the splitting idea |35|, |36(. We distinguish 
incoming and outgoing waves along the boundaries for the linear kinetic subproblem, and 
then unite the potential energy subproblem to yield nonlinear boundary conditions which 
are of local form and examine the coupling between the equation governing the wave in 
the computational domain and the boundary conditions on the boundary yields a well- 
posed problem. Furthermore, the obtained boundary conditions are "easy" to discretize and 
"cheap" to compute in terms of the computational time. We will then perform numerical 
tests to illustrate the effectiveness and tractability of this approach. 

II. NONLINEAR ABSORBING BOUNDARY CONDITIONS 
A. Linear Schrodinger equation 

Let us first consider the local ABCs of the linear Schrodinger equation, 

ih W(x,t) = (2) 
dt 2m dx 2 ' 

which will be the basis of constructing ABCs of nonlinear problems, because it is the kinetic 
part of the NLS equation ([!]). 
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In the frequency domain of the Fourier transform, the one way equations at the two 
boundaries, which annihilate all the outgoing waves, can be represented by 



i\j 7^?Ar ± v 7 ^ = 0, (3) 

where the subscript represents the partial derivative, ip(x,u) is the Fourier transform of 
i(j(x, t), oj is the frequency, the plus sign in "±" corresponds to the right boundary condition, 
and the minus sign corresponds to the left. The transparent boundary conditions in exact 
manner of the problem is then derived by an inverse Fourier transform: 



, h , e 4 ! d f* , . , dr 

*Vs^ ± W*y.* (l,T) ^7 = 0, (4) 

which is a nonlocal boundary condition. In order to get local boundary conditions, a usual 
way is to approximate the square root -Jui by using polynomials or rational polynomials. 
For instance, applying the approximations 

u = ^/u^+0(u-Uq), (5) 



ZJ=^ + -^ + O(c-c ) 2 , (6) 
- v ^(3o; + a;o) 3 

u = * — + 0(w - uj ) , 7 

uj + 3u> 

we obtain the first three absorbing boundary conditions in the physical domain: 



1st order : i\j — ± k^ip = 0, 

2m ox 

2h, dip , , , 




2nd order : iip t ± i\ — ko— — h k Q ip = 0, 

V m ox 

border: - ± + ^ = . (8) 



m dx dxdt \ m dt J 

Here, uj is a positive constant, and k = ^JoJ^ which is called the wavenumber parameter. 



B. Nonlinear absorbing boundary conditions 

In this section, we derive the absorbing boundary conditions for the one-dimensional NLS 
equation (JTJ) on the region Q = [xi,x r ] by assuming that there do not exist outgoing waves 
at the two ends X\ and x r . We remark that this assumption is not correct for a general 



wave due to the nonlinear or space-dependent potential 



40j, but it is reasonable for many 



situations of physical applications; in fact, this is a fundamental assumption so that we can 
design absorbing boundary conditions. 

To understand the philosophy of obtaining nonlinear absorbing boundary conditions, let 
us rewrite Eq. ([T]) into the following operator form: 

•R^^ = [f + Vty(*,t), (9) 

where T represents the linear differential operator which accounts for the kinetic energy 
part, and V is the nonlinear operator that governs the effects of the potential energy and 
nonlinearity. These are given by 

and V = V(x) + f(\^). (10) 

In a time interval from t to t+r for small r, the exact solution of Eq. (jUJ) can be approximated 
by 

t + r) « e- i[f+ ^ T/n ^(x, t) « e" ifT/ V^ T// V(x, *). (11) 

This means in a small time step r that the approximation carries out the wave propagation 
in a kinetic energy step and a potential energy step separately. This is the basic idea of the 
well-known time-splitting method (or the split-step method) which is effective to numerically 
solve the nonlinear Schrodinger-type equations. 

Since we have assumed that there are no outgoing waves on the artificial boundaries 
xi and x r , noticing formulae ©, the kinetic operator T can be approached by one-way 
operators 

^ ,2, - fi (Wih§ + 4 (12) 

t « f w = -tf ( 2i|- ± ak ) "75iM|_ ± My (13) 

\ ox J \ m ox m J 
which correspond to the second- and third- order local boundary conditions in Eq. (jSj), 
respectively. Here, we also take the positive sign in "±" at x = x r and take the negative 
sign at X[. Substituting T^ n ' for T in Eq. f fTTl) yields approximating operators to the solution 
operator as follows, 

which imply the one-way equations, 

ih^^ = [f^+V]^(x,t). (15) 



They approximate to Eq. (jUJ) at two boundaries and being going to act as absorbing bound- 
ary conditions we need. Concretely, we obtain nonlinear absorbing boundary conditions: 



n 



n 



ih 



dip 
dt 



I2h 8 

< ±l S^Wx +kl) + V{x) + fm2) 



(16) 







ih ( 2i— — ± 6k 

ox 



dip 
~dt 



ft 2 f f) 

dx 



m 



+{V(x)+f(\iP\ 2 )}(2i—±6k 







ip. 



(17) 



Here, in order to increase the accuracy of the boundary conditions, the wavenumber param- 
eter kn is a function of time t, which will be determined by a windowed Fourier transform 
I — I 

method [36] based on the frequency property of the wave function near the artificial bound- 
aries. We express the boundary conditions f)16p and (j!7p in the operator forms: 

M { "\x,t)ip(x,t) = 0, (18) 
M { +\x,t)ip{x,t) = 0, (19) 

where B_ and B + respectively represent the left and right boundary conditions; that is, the 
minus sign in " ± " is taken for B_, and the plus sign for B + . The initial value problem on 
the open domain of Eq. (CD) restricted to the truncated interval [ then approximated 

by an initial boundary value problem with local absorbing boundary conditions: 



dip(x, t) 
~dt 



h 2 d 2 



2m0x 2+V ^ + f ^ 
ip(x,0) = ipo(x), for x&[xi,x r ], 

M^\x u t)il)(x u t) = 0, for n= 2, or 3, 



ip, for x G [xi, x r 



(x r , t)ip(x r , t) = 0, for n = 2, or 3. 



(20) 
(21) 
(22) 
(23) 



C. Well-posedness of the induced IB VP 

It is natural to examine the well-posedness that we concern most of the induced ini- 
tial boundary value problem confined in the finite computational interval for the nonlinear 
Schrodiger equation by using the nonlocal absorbing boundary condition. The classical 
energy method seems to be difficult to obtain the stability result of the problem under 
consideration here. In this section, we will investigate the stability based on Kreiss's nor- 

26J and references therein, which was also used for the wave 



mal mode analysis method [14], 
equations and the linear Schrodinger equation. 



The idea of the general algebraic normal mode analysis is based on the fact that the well- 
posed problem does not admit any complex eigenvalues with positive real parts Re(s) > 0, 
or generalized eigenvalues with Re(s) = and the negative (positive) group velocity of 
normal mode on the right-hand (left-hand) boundary. The eigenvalues and generalized 
eigenvalues satisfy the dispersion relations of both the interior equation and the equation 
on the boundaries. If there exist such eigenvalues, the problem will admit an unboundedly 
grown normal mode e st and is hence unstable. If there exit such generalized eigenvalues, 
the boundary conditions will admit a incoming wave which will propagate energy into the 
interior domain to disrupt the solution in the computational interval and generate a spurious 
wave solution. 

We need to check if there exist eigenvalues and generalized eigenvalues by replacing the 
solution of the plane wave form ip(x,t) = e st+kx into Eq. (120]) and boundary conditions 
(I22p (I23p . There are left-hand and right-hand boundary, they can be considered by the same 
argument. Here we only consider the wellposedness with the right-hand boundary conditions. 
For n=2, substituting the plane wave into the dispersion relation and the boundary condition 
yields 



In very small time step r the basic idea of the well-known time-splitting method is to carry 
out the wave propagation in a kinetic energy step and a potential energy step separately. 
When time step is taken limit to zero, that is the reason why we can write the nonlinear 
shrodinger equation in above forms, the nonlinear term is considered as a potential function. 
Hence, let (J23J minus f[2"5l) and substitute the result into (T21]) to arrive at 



where V + f — V\ + fi + i(v 2 + f 2 )- For the 3rd order approximation on the right boundary, 
by the same argument, we have 




(24) 



(25) 




(26) 




(28) 



(27) 
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Combine (1271) and (1251). we have 



2m 2m ik + 3&o 
solve (1291) and substitute the results k = ik into (127)) to obtain 

. = ^_WEL£M (30, 

•=- 

It is easy to have Re(s) < for v 2 + f% < from ( 126)) and ( 13T1) whatever n = 2 or n = 3. 

In particular, when 1/ and / are both real functions, i.e. U2 + /2 — 0, we thus have Re(s) = 
0, hence s is wholly imaginary. Since the boundary condition is designed to annihilate all the 
outgoing wave, this implies that the group velocity is positive on the right-hand boundary. 
Hence there is no generalized eigenvalue which will propagate energy into the interior interval 
to influence on the true wave solution. If any instabilities which might be admitted by the 
generalized eigenvalues of the ABCs exit, they would not propagate wave to affect the interior 
solution. Keep this concept in mind, by linearizing the nonlinear terms in a sufficiently small 
time step, we can say that the boundary conditions of the problem ( 1201) - (1231) are well-posed 
if functions V(-) and /(•) satisfy v 2 + f% < 0. Hence the boundary condition for n=2 or 3 is 
well-posed and numerical tests below show that the boundary condition is stable. 

III. NUMERICAL EXAMPLES 

In this section we will test the efficiency of the nonlinear absorbing boundary conditions 
by using the finite-difference scheme to the induced initial boundary value problem (120)) - 
023]). 

In the computational interval [xi, x r ], let J is a positive integer, denote Ax = (x r — x-j)/ J and 
At by grid sizes in ^space and time, respectively. We use the following Crank-Nicolson-type 
difference scheme 



11 sum 



lH At = -2^ A^ + ^ + m 2 1 ^ 



T 2 
3 ' 



(32) 



to discretize the nonlinear Schrodinger equation inside the computational domain, which is 
unconditionally stable and of second order in time. Here ifjj 2 = — — ^ — - and ip™ represents 

8 



the approximation of wave function ip at the grid point (xj,t n ) with j = 0,1,..., J, xq = 
xi, Xj = xi + jAx, xj = x r , and t n = nAt. Obviously, the above systems of equations can 
not be solved uniquely since the number of equations are less than that of unknowns. At 
the two unknown ghost points X-\ and xj + \ no equation can be defined from (J22D • Thus 
the absorbing boundary conditions are required to provide two extra equations on the grid 
points x s for s G S = { — 1, 0, 1} U {J — 1, J, J + 1} in the vicinity of boundary. The 3rd 
order ABCs are discretized as follows: 



°£&-2i(V a + f(\ti + t\*) 



2 Ax Ax At 1 Ax At 



±6M^^ ± - QhV s - 6fc /(|V>r¥)) ^ = 0, (33) 

where the positive sign in " ± " corresponds to the right boundary s = J and the negative 

one corresponds to the left boundary s = 0. The above scheme is implicit, we need use an 

n , i 

iterative strategy to solve the nonlinear equations (I32l - fl33l by replacing f(\ipj 2 | 2 ) with 
f(\~ — 2 ^l 2 )' w here the superscript k denotes the k-ih iteration at each time step and 



j 



the initial iteration is given as ("ip™ +1 )° = ip 
Without loss of generality, we rewrite the original NLS equation (JTJ) as 

#t = -^xx + ^/(|^| 2 )^ + ^, (34) 

where gf{\ip\ 2 )ip represents different physical phenomena when chosen different nonlinear 
forms. In example 1, we take f{\ip\ 2 ) = \ip\ 2 , V = and g = —2, which corresponds 
to a focusing (g < 0) effect of the cubic nonlinearity. Many applications in science and 
technology have connections with this kind of NLS equations |27|, |30j. Furthermore, in 
this case equation (|34|) can be solved analytically by the inverse scattering theory, its exact 
soliton solution has the following form 



if)(x, t) = A^-^sech(Ax - 2ABt)e lBx+i{A2 - B2) \ (35) 

where A and B are real parameters:^ determines the amplitude of the wavefield and B is 
related to velocity of the soliton. In example 2, we take fi\ip\ 2 ) = |^| 4 , V = and g = —2, 
the NLS equation becomes of a quintic nonlinearity. The computational interval is set to be 
[xi,x r ] = [—5,5]. For this kind of NLS equation with power-like nonlinearity, if the initial 
energy satisfies: 

E(if> ) := \\VM\h-l\\M%<0, (36) 



TABLE I: The values and orders of Li-norm at different times and mesh sizes. 





Ax = 0.2 


order 


Ax = 0.1 


order 


Ax = 0.05 


order 


t = 2 


1.570e-2 




4.149e-3 


2.076 


1.017e-3 


2.029 


t = 3 


2.946e-2 




6.797e-3 


2.116 


1.656e-3 


2.037 


t = 4 


3.622e-2 




8.091e-3 


2.162 


1.955e-3 


2.049 


t = 5 


3.145e-2 




6.948e-3 


2.178 


1.724e-3 


2.011 


t = 6 


2.691e-2 




6.010e-3 


2.163 


1.581e-3 


1.926 



the solution ip will blow up in finite time ([5j] and references therein), i.e, there exists T > 
such that 

lim 1 1 V x ib\ I ip. — > +oo. 

t->T 

In example 3, we take f{\ip\ 2 ) = \ip\ 2 , V ^ and g = 2, the NLS equation represents a 
nonlinear wave with repulsive interaction. 

Example 1. We first consider the case of g = —2 and the temporal evolution of an initial 
bright solition 

ip(x, 0) = sech(x - x )e 2i{x - Xo) (37) 



with xo = 15. The exact solution is the form of (|35|) with A = 1 and B = 2. This case is 
used to test the performance of the nonlinear ABCs when the wave impinges on the right- 
hand-side boundary of the computational interval [0,30]. FIG. 1 shows the evolution of 
absolute error \ip e xa — VVil on ^ ne artificial boundary point by picking wavenumber ko = 2, 
spatial size Ax = 0.1 and time size At = 0.01 and its refinement mesh. When refining the 
mesh, one can see that the error will converge quickly. In order to explore the convergence 
rates in space and time, we take time size as At = Ax 2 . Denote the L\— norm by 

1 J N 

El = (JTWnTT) £ £ W ~ nv (38) 

\ ~ )\ ~ > j=o n= Q 

FIG. 2 shows the evolution of the L\— norm with different mesh sizes. From this figure, we 
observe a second-order accuracy of the errors. To character convergence rates exactly, we 
take the different values of the Li-norm at different times and mesh sizes. TABLE I shows 
the second-order convergence rate. 

Another im po rt a„ t m etho d Q to 8 ee the infl ue n ee of t a on the bo UIldary conditions is 



10 



TABLE II: The values and orders of Li-norm at different times and mesh sizes. 





0.5 


0.75 


1.0 


1.25 


1.5 


1.75 


2.0 


2.125 


2.25 


-2 


6.07e-2 


1.601e-2 


4.70e-3 


1.52e-3 


5.17e-4 


1.81e-4 


7.30e-5 


6.13e-5 


7.06e-5 


-10 


6.07e-2 


1.60e-2 


4.70e-3 


1.52e-3 


5.17e-4 


1.81e-4 


7.30e-5 


6.13e-5 


7.05e-5 


g\k 


2.5 


2.75 


3.0 


3.25 


3.5 


3.75 


4.0 


4.25 


4.5 


-2 


1.51e-4 


3.23e-4 


6.10e-4 


1.04e-3 


1.62e-3 


2.40e-3 


3.39e-3 


4.60e-3 


6.06e-3 


-10 


1.51e-4 


3.23e-4 


6.10e-4 


1.03e-3 


1.62e-3 


2.40e-3 


3.39e-3 


4.60e-3 


6.06e-3 


g\h 


4.75 


5.0 


5.25 


5.5 


5.75 


6.0 


6.25 


6.5 


6.75 


-2 


7.77e-3 


9.73e-3 


1.20e-2 


1.44e-2 


1.72e-2 


2.02e-2 


2.34e-2 


2.69e-2 


3.05e-2 


-10 


7.77e-3 


9.73e-3 


1.20e-2 


1.44e-2 


1.72e-2 


2.02e-2 


2.34e-2 


2.69e-2 


3.05e-2 



the reflection ratios. The reflection ratios [25j are calculated as follow 



(39) 



which is efficient to measure the performance of nonlinear ABCs. r = implies that the 
soliton passes through the boundary completely, and r = 1 indicates that the wave is 
completely reflected into the computational domain. In the experiments of BECs, g often 
takes larger value so that the nonlinear term will become very strong. TABLE II presents 
the reflection ratios with parameters g being chosen -2 and -10, respectively. We observed 
the wave numbers influence the reflection ratios under spatial size Ax = 0.1 and time 
size At = 0.05. But the reflection rates r are insensitive to ko, which can be taken in a 
large interval [1.0,5.0] such that the reflection ratios is less than 1%. Out of the interval 
r will increase rapidly. Furthermore, whatever the parameter g is given large or small, i.e. 
whatever the nonlinear term is strong or not, the reflection ratios are insensitive to the 
nonlinear term, either. The parameter g can even be chosen in a larger interval. FIG. 3 
plots the amplitude of the soliton under ko = 2.0 and spatial size Ax = 0.1, respectively. 
No reservable reflections can be detected at all. 

Example 2. We consider the quintic NLS equation with the initial value 

where the wavenumber ko = 8. A direct calculation obtains the initial energy E(if)o) « 
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80.5478 > 0, which implies that the numerical experiment will not blow-up. By using the 
Crank-Nicolson scheme proposed by Duran et al 7], FIG. 4 shows the reflection ratios for 
different wavenumbers under spatial size Ax = 0.1 and its refinement mesh size Ax = 0.05, 
respectively. Comparing the two curves, we can see that the reflection ratios decrease fast 
when wavenumber taken in the vicinity of the initial wavenumber ko = 8. And the reflection 
ratios in larger mesh is larger than these in smaller mesh when ko < 4. This indicates that 
the parameter ko in nonlinear ABCs should be given in suitable way, i.e, we had better endow 
the wavenumber with being equal to the initial wavenumber if known. We can also see if 
the ko is equal to half of the group velocity of the wave impinged on the boundary point for 
third-order case, the nonlinear ABCs will not disrupt the the interior solution when mesh 
is refined and refined. FIG. 5 shows the wave amplitude of utilizing the Crank-Nicolson 
method proposed by Duran under mesh size Ax = 0.01, At = 0.001 and ko = 8.0. The 
reflecting wave can not be observed and the reflection ration is 7.637e-5. 
Example 3. We consider a nonlinear wave with repulsive interaction with g = 2 in ( l34j ). 
The initial data and potential function are taken to be Gaussian pulses 

i/j(x, 0) = e -°- l{x - xo)2 and V(x) = e~°- 5{x - xo)2 (40) 



with xq = 15. It has been an example in 35|, |36), and is used to model expansion of a Bose 



Einstain condensate composed of waves with different group velocities. In the calculation, 
we choose L = 30, Ax = 0.1 and At = 0.0375. FIG. 6 (a)-(b) depict the evolution of the 
wave's amplitude with the third-order nonlinear ABCs for different wavenumbers ko. We 
emphasize the very small spurious reflections by highlighting the figures, which can make 
us visualize them by showing the associated shadow zones. And the reflective waves are too 
small to be seen with wavenumbers chosen in the appropriate range [1.25,2]. In FIG. 6 (a) 
and 6(b), no observable reflection can be defected at all. If we take the boundary conditions 



as Dirichlet or Neumann boundary condition, the reflective wave is very large (ref. to 351 ]) 
at t=6 with the same meshes. 

We remark that picking the suitable wave-number k in the nonlinear ABCs plays a very 
important role if we want to obtain an efficient nonlinear absorbing boundary condition 
with more and more refinement mesh. The ideal wavenumber ko is equal to a half of group 
velocity when a single solitary wave impinges on the boundary. For an arbitrary wave 



packet, a suitable wavenumber can be picked adaptively by using the Gabor transform 
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36] 



to capture local frequency information in the vicinity of artificial boundary, for which we 
will not discuss in detail in this paper. 

IV. CONCLUSION 

Through approximating the kinetic energy term by a one-way equation and uniting it 
with the potential energy equation, we proposed an efficient method for designing absorbing 
boundary conditions of nonlinear Schrodinger equations. The obtained boundary condi- 
tions, which are of nonlinear form, can perfectly absorb the outgoing waves. We examined 
the coupling between the equation governing the wave in the computational domain and 
the boundary conditions on the boundary yielded a well-posed problem in a sufficiently 
small time interval. Furthermore, some numerical examples are given to demonstrate the 
effectiveness and tractability of the artificial boundary conditions. 
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FIG. 3: Evolution of the amplitude of the bright soliton with ko = 2.0. 




FIG. 4: The reflection ratios when chosen different wavenumbers ko under mesh Ax = 0.1 and its 
refinement mesh. 




FIG. 5: Evolution of the amplitude of the wave for the quintic NLS equation. 
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